Table 1 (conV) 



* printO ~ only routine visible outside this module 

5 * static: 

* getmatO - trace back best path, count matches: printQ 

* pr_alignO - print alignment of described in array pQ: printO 

* dumpblockO ~ dump a block of lines with numbers, stars: pr_align0 

* numsO ~ put out a number line: dumpbiockO 

10 * putlineO - put out a line (name, [num], seq, [numl): dumpblockO 

* starsO - -put a line of stars: dumpblockO 

* stripnameO - strip any path and prefix from a seqname 
*/ 



Mefine SPC 3 

#defineP_LINE 256 /* maximum output line */ 

#define P SPC 3 /* space between name or num and seq */ 

extern _day[26][26]; 

int olen; /* set output line length */ 

FILE *&; /* output file */ 

25 printO print 
{ 

int Ix, ly, firstgap, lastgap; /* overlap */ 

if ((fe = fopen(ofile, "w")) = = 0) { 
30 ^rintf(stderr,"%s: can't write %s\n", prog, ofiie); 

cleanup(l); 

} 

ftirintf(fs, "< first sequence: %s Gength = %d)\n", namex[0]. lenO); 
^rintf(fx, "< second sequence: %s (length = %d)\n", namex[l], lenl); 
35 olen = 60; 

Ix = lenO; 
ly = lenl; 

firstgap = lastgap = 0; 

if (dmax < lenl - 1) { /* leading gap inx */ 
40 pp[0].spc = firstgap = lenl - dmax - 1; 

ly-=pp[0].spc; 

> 

else if (dmax > lenl - 1) { /* leading gap in y */ 
pp[l].spc = firstgap = dmax - (lenl - 1); 
45 Ix-=pp[l].spc; 
} 

if (dmaxO < lenO - 1) { /* traUing gap in x */ 
lastgap = lenO - dmaxO -1; 
Ix -= lastgap; 

50 } 

else if (dmaxO > lenO - 1) { /* traili^ gap in y */ 
lastgap = dmaxO - (lenO - 1); 
ly -= lastgap; 

55 getmatOx, Ij 

pr_alignO; 
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Table 1 (conV) 

I* 

* trace back the best path, count matches 

*/ 

static 

5 getmat(lx, ly, firstgap, lastgap) getmat 
int Ix, ly; /* "core" (minus endgaps) */ 

int firstgap, lastgap; /* leading trailing overlap */ 

int nm, iO, il, sizO, sizl; 

10 char outx[32]; 

double pet; 

register nO, nl; 

register char *pO, *pl; 

15 /* get total matches, score 

*/ 

iO = il = sizO = sizl = 0; 
pO = seqx[0] + pp[l].spc; 
pi = seqx[l] + pp[0].spc; 
20 nO = pp[l].spc + 1; 

nl = pp[0].spc + 1; 

nm = 0; 

while (*pO&&*pl){ 
25 if (SizO) { 

pi++; 

nl+ + ; 
sizO-; 

30 else if (sizl) { 

pO++; 



nD++; 
sizl—: 



if (xbm[*pO-'A']&xbm[*pl-'A']) 

nm++; 
if (nO++ ==pp[0].x[iO]) 

sizO = pp[0].n[iO++]; 
if(nl + + ==pp[l].x[il]) 

sizl = pp[l].n[il + +]; 

pO+ + ; 
pl + +; 



/* pet homology: 

* if penalizing endgaps, base is the shorter seq 

* else, Imock off overhangs and take shorter core 

*/ 

if (endgaps) 

Ix = (lenO < lenl)? lenO : lenl; 

else 

Ix = Ox < ly)? Ix : ly; 
pet = 100.*(double)mii/(double)lx; 
^rintf(fe, "\n"); 

Q)rintf(fx, " < %d matches in an overlap of %d: %.2f percent similarityVn", 
nm, (nm == 1)? : "es", Ix, pet); 
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' Table 1 (conf) 

^rintf(fx, " <gaps in first sequence: %d", gapx); getmat 
if(gapx){ 

(void) spriDtf(outx, " (%d %s%s)", 

i^apx, (dna)? "base": "residue", (ngapx == 1)? "":"s"); 

^rititf(fe,"%s", outx); 

^rintf(fe, gaps in second sequence: %d", gapy); 
if(gapy){ 

(void) sprintf(outx, " (%d %s%s)", 

ngapy, (dna)? "base":"residue", (ngapy == 1)? "":"s''); 

lprintf(fic,"%s", outx); 

} 

if (dna) 

^rmtf(fx, 

"\n< score: %d (match = %d, mismatch = %d, gap penalty = %d + %d per base)\n", - 
smax, DMAT, DMIS, DINSO. DINSl); 

else 

Qirititfiffe, 

"\n< score: %d (Dayhoff PAM 250 matrix, gap penalty = %d + %d per residue)\n", 
smax, PBSrSO, PINSl); 
if (endgaps) 

^rintf(fx, 

"<en^aps penalized, leftendgap: %d %s%s, right endgap: %d %s%s\n", 
firstgap, (dna)? "ba^" : "residue", (firstgap == 1)? "" : °s", 
lastgap, (dna)? "base" : "residue", (lastgap == l)? : "s"); 



else 



Q>rintf(fic, "< endgaps not penalized\n"); 





stadc 






stafic 


hnax; 




static 


ij[2]; 


35 


static 


nc[2I; 


static 


ni[2]; 




static 
static char 


siz[2]; 




static char 


*ps[21; 
*po[2]; 


40 


static char 


out[2][P_LINE]; 


static char 


star|P_LINE]; 



/* matches in core ~ for checking */ 
/* lengths of stripped file names */ 
/* jmp index for a path */ 
/* number at start of current line */ 
/* current elem number — for gapping */ 

/* ptr to current element */ 
/* ptr to next output char slot */ 
/* output line */ 
/* set by starsO */ 

/* 

* print alignment of described in struct path ppO 
static 

pr_alignO pr_aligll 
tot nn; /* char coimt */ 



45 



for (i = 0, Imax = 0; 1 < 2; i+ +) { 
nn = stripname(namex[i]); 
if (nn > Imax) 
55 Imax = nn; 

nc[i] = 1; 
nipl = 1: 
siz[i] = ij[i] = 0: 
OO ps[i] = seqx[i]: 

pop] = out[i]; 
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